%% 做出一个随机旋转角度的连接三角

% function [输出变量] = 函数名称(输入变量）
function [x2,y2,locm,triangle1,triangle2] = creatanconnectedTriangle(triangle_radius,x1,y1) %x1,y1是第一个椭圆的中心点位置
% x2，y2第二个矩形中点坐标（画线）；locm刺激中点位置（控制刺激距离）；rect矩形坐标

[xax, yax] = meshgrid(1:227, 1:227);
rax = xax+1i*yax;

loc1 = x1+1i*y1;

% 三角2中心位置
rt1 = abs(rax-loc1);%rt1<circle_radius就是圆1;下行命令是rt1<4*r，就是为了找出圆2运动轨迹围成的大圆，然后再用"bwboundaries"抠出轨迹
loc2ind = rt1<4*triangle_radius; % 第2个圆的圆心运动轨迹的大圆
loc2bound = bwboundaries(loc2ind); % 第2个圆的圆心运动的轨迹（loc2ind的圆周长），是cell
loc2bound = loc2bound{1};% 将cell提取出来，形成矩阵，方便下面提取
loc2 = loc2bound(randi(size(loc2bound,1),1),1:2); % 提取圆2的圆心运动坐标的范围，随机挑一个坐标作为圆2的坐标
loc2 = loc2(2) + 1i*loc2(1);% 实部对应x（列），虚部对应y（行）

% 三角2的中点坐标参数
x2 = real(loc2); % 圆心2的横坐标
y2 = imag(loc2); % 圆心2的纵坐标

% 刺激中点位置
xm = (x1+x2)/2;
ym = (y1+y2)/2;
locm = xm + 1i*ym;

% 三角坐标位置
% 第一个三角
ux1 = x1; % 上顶点x坐标
uy1 = y1-triangle_radius; % 上顶点y坐标
lx1 = x1-sqrt(3)/2*triangle_radius; % 左顶点x坐标
ly1 = y1+1/2*triangle_radius; % 左顶点y坐标
rx1 = x1+sqrt(3)/2*triangle_radius; % 右顶点x坐标
ry1 = y1+0.5*triangle_radius; % 右顶点y坐标
triangle1 = [ux1,uy1;lx1,ly1;rx1,ry1];

%第二个三角
ux2 = x2; % 上顶点x坐标
uy2 = y2-triangle_radius; % 上顶点y坐标
lx2 = x2-sqrt(3)/2*triangle_radius; % 左顶点x坐标
ly2 = y2+1/2*triangle_radius; % 左顶点y坐标
rx2 = x2+sqrt(3)/2*triangle_radius; % 右顶点x坐标
ry2 = y2+0.5*triangle_radius; % 右顶点y坐标
triangle2 = [ux2,uy2;lx2,ly2;rx2,ry2];


end